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Abstract 

The Laplacian growth problem in the limit of zero surface tension is proved 
to be equivalent to finding a particular solution to the dispersionless Toda lattice 
hierarchy. The hierarchical times are harmonic moments of the growing domain. 
\ The Laplacian growth equation itself is the quasiclassical version of the string 

equation that selects the solution to the hierarchy. 
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The Laplacian growth problem is one of the central problems in the theory of pattern 
formation. It has many different faces and a lot of important applications. In general 
words, this is about dynamics of moving front (interface) between two different phases. In 
many cases the dynamics is governed by a scalar field that obeys the Laplace equation; 
that is why this class of growth problems is called Laplacian. Here we shall confine 
ourselves to the two-dimensional (2D) case only. To be definite, we shall speak about 
two incompressible fluids with different viscosities on the plane. In practice, the 2D 
geometry is realized in the narrow gap between two plates. In this version, this is known 
as the Saffman- Taylor problem or viscous fingering in the Hele-Shaw cell. For a review, 
seei. 

We shall mostly concentrate on the external radial problem for it turns out to be 
the simplest case in the frame of the suggested approach. Let the exterior of a simply 
connected domain on the plane be occupied by a viscous fluid (oil) while the interior 
be occupied by a fluid with small viscosity (water). The oil/water interface is assumed 
to be a simple analytic curve. Other versions such as internal radial problem, wedge or 
channel geometry are briefly discussed at the end of the paper. Basically, they allow for 
the same approach. 

Let p(x,y) be the pressure, then p is constant in the water domain. We set it equal 
to zero. In the case of zero surface tension p is a continuous function across the interface, 
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so p = on the interface. In the oil domain the gradient of p is proportional to local 
velocity V = (V x , V y ) of the fluid (Darcy's law): 

V = — Kgradp, (1) 

where k is called the filtration coefficient^. In particular, this law holds on the interface 
thus governing its dynamics: 

V. = -«|. (2) 

Here V n is the component of the velocity normal to the interface and dp/dn is normal 
derivative. Since the fluid is incompressible (divy = 0), the Darcy law implies that the 
potential <&(x,y) = —np(x,y) is a harmonic function in the exterior (oil) domain: 

A$(x,|/) = 0, (3) 

where A = d% + dy is the Laplace operator on the plane. The asymptotic behaviour of 
the function $ very far away from the interface (at infinity) is determined by the physical 
condition that there is a sink with constant capacity q placed at infinity. This means 
that 

(V, n)dl = q . 



where 7 is any closed contour encircling the water domain, n is the unit vector normal 
to 7. So, we require $ = on the interface and $ = ^Tog \z\ as \z\ — > 00. The goal is to 
describe the interface motion subject to the local dynamical law V n = dQ/dn. 

An effective tool for dealing with this problem is the time dependent conformal map- 
ping technique (see e.g. [Q). Passing to the complex coordinates z = x + iy,z = x — iy 
on the physical plane, we bring into play a conformal map from a reference domain on 
the mathematical plane w to the growing domain on the physical plane. By the Riemann 
mapping theorem, such a map does exist and, under some conditions, is unique. More 
precisely, let z = \(w) be the univalent conformal map from the exterior of the unit 
circle to the exterior of the interface (i.e., to the oil domain) such that 00 is mapped 
to 00 and the derivative A'(oo) is a positive real number r. Under these conditions the 
map is known to be unique. The Laurent expansion of the X(w) around 00 has then the 
following general form: 

00 

X(w) = rw + UjW~ 3 . (4) 

If the interface moves, the conformal map z(w,t) = \(w,t) becomes time-dependent. 
The interface itself is the image of the unit circle \w\ = 1: as w = e 1 ^, < <p < 2tt, sweeps 
over the unit circle, z = X(e l< ^,t) sweeps over the interface at the moment t. 

Having defined the conformal map, we immediately see that the real part of the loga- 
rithm of the inverse conformal map, w(z), provides the solution to the Laplace equation in 
the oil domain with the required asymptotics: $(x,y) = ^Relogw(z). Let us introduce 
the complex velocity V = V x — iV y . The obvious formula 

q d\ogw(z,t) 9$ 9$ 

%- 



2tt dz dx dy 



1 This coefficient is inverseiy proportional to the viscosity, so the Darcy law is formally valid in the 
water, too. 
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allows one to represent the Darcy law ([I]) as follows: 

q dlogw 

v{z) = ~^z~ ' (5) 

where the derivative is taken at constant t. 

In terms of the time-dependent conformal map, the Darcy law is equivalent to the 
following relation referred to as the Laplacian growth equation (LGE): 

(dz dz \ q 

It first appeared in 1945 [@] in the works on the mathematical theory of oil production. 
^From now on we set q = tt without loss of generality. (This amounts to a proper rescaling 
of t.) Introducing the Poisson bracket notation 

r f df dg dg df 

{f > 9} = W d^dt- W d^dt> (7) 
for functions / = f(w, t), g = g(w, t) of w, t, we rewrite the LGE in the suggestive formF] 

{z(w,t),z(w-\t)} = l. (8) 

The LGE thus means that the transformation from log to, t to z, z is canonical. 

For a technical simplicity, we assume that the point z = lies in the water domain. 
The interface dynamics given by the Darcy law (or, equivalently, by the LGE) implies 
that if a point (x, y) is in the water (interior) domain at the initial moment, then it 
remains there for all values of time. In particular, our assumption that the point z = 
belongs to the water domain means that zeros of the function A (it;) are inside the unit 
circle for any t. 

The Laplacian growth is a particular case of the 2D inverse potential problem. The 
shape of the interface can be characterized by the harmonic moments Ck of the oil domain 
and the area Co of the water domain: 

C k = - 1 1 z~ h dxdy, k>l; C = [ f dxdy . (9) 

J ./exterior J ./interior 

(The integrals at k = 1, 2 are assumed to be properly regularized.) A remarkable result of 
Richardson || shows that the LGE (f|) implies conservation of the harmonic moments Ck 
when the interface moves, dC^/dt = 0, while area of the water domain grows linearly in 
time: Co = id. Therefore, the problem can be posed as follows: to find the shape of the 
domain as a function of its area provided all the harmonic moments of the exterior are 
kept fixed. Since the harmonic moments are coefficients in the expansion of the Coulomb 
potential created by a homogeneously distributed charge in the oil domain, this is just a 
specification of the inverse potential problem. We remind that to know the shape of the 
domain is the same as to know the coefficients r, Uj of the conformal map (f|). 

A good deal of hints that the LGE has much to do with integrable systems have 
been known for quite a long time @-|J. However, its status in the realm of integrability 

2 Given a Laurent series f{z) — ^2jfjZ J , we set f(z) = ^2jfjZ J > so z(w) and are complex 

conjugate only if |w| = 1. 



3 



was obscure until very recently. The nature of the LGE and its relation to integrability 
are clarified in the work 0. The idea is to treat the LGE not as a dynamical equation 
but as a constraint in a bigger integrable hierarchy. The latter turns out to be an 
infinite Whitham hierarchy of the type first introduced in ||. This hierarchy is a multi- 
dimensional generalization of integrable hierarchies of hydrodinamic type JJ. It naturally 
incorporates the general inverse potential problem as well. Namely, the coefficients of 
the conformal map as functions of all the harmonic moments are given by a particular 



solution to the dispersionless Toda lattice hierarchy (see |L0[ for a detailed study of 
the latter). The hierarchical evolution times are just the harmonic moments and their 
complex conjugate. In the Laplacian growth all of them but the area Cq are frozen. 
Making them alive, one moves over the space of initial data for the LGE, and recovers 
the Whitham hierarchy. 

For a more convenient formulation of the result, let us rescale the harmonic moments 
and introduce the new notation for them: 

t = t = — , t k = % i k = % k>l. (10) 

7T TXK TTK 

The symbols {f{w))± below mean a truncated Laurent series, where only terms with 
positive (negative) powers of w are kept, (f(w)) is a constant part (w°) of the series. 

Theorem 1 The conformal map (fjj obeys the following differential equations with re- 
spect to the harmonic moments: 

^)= W ,A(»-)}, *5 = - { 4*-)}, (12) 
where the Poisson bracket is defined in H), and 

H j (w) = (X(w)) + + ±(X(w)) o , (13) 
H j (w)=(\i(w- 1 ))_ + ±(\i(w- 1 )) o . (14) 



The proof is sketched in |7]] and [|Tl|]. These are the Lax-Sato equations for the dispersion- 
less Toda lattice hierarchy of non-linear differential equations, with the X(w) and X(w _1 ) 
being the Lax functions. On comparing coefficients in front of powers of w in ([11]), (03), 
one obtains an infinite set of non-linear differential equations for the coefficients of the 
comformal map. Altogether, they form the hierarchy. The particular solution that solves 
the inverse potential problem is selected by the constraint (§), where z{w,t) = X(w,t), 
z(w, t) = A(u> -1 , t), and all tk are fixed. This constraint is known as (a quasiclassical ver- 
sion of) the string equation. Surprisingly, this very constraint is the key ingredient of the 



integrable structures in 2D gravity coupled with c = 1 matter [|12| . |T3f . The mathematical 
theory of the dispersionless hierarchies constrained by string equations was developed in 
]b|, [L5| and extended to the Toda case in [|Kj . 
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The Lax-Sato equations imply |1JJ the existence of the prepotential function F(t, t k , i k ). 
This function solves the inverse potential problem in the following sense. Let C_ k , k > 1, 
be the complimentary set of harmonic moments, i.e., the moments of the interior of the 
domain: 

C_ k = [f z k dxdy. (15) 
J ./interior 

Were we able to find from a given set C , C k (i.e., given t, t k ), this would yield 
the complete solution, alternative but equivalent to knowing the conformal map @. 
Indeed, using the contour integral representation of the harmonic moments, it is easy to 
see that the generating function of all the harmonic moments, obtained as an analytic 
continuation of the Laurent series 

1 OO 1 oo 

//(a) = - e c k \ k = E kt ^ k + 1 + - E c -^~ k , (is) 

n fcez k=l 71 k=l 

would allow us to restore the interface curve via the equation \z\ 2 = fi(z) (S(z) = \i{z)jz 
is what is called Schwarz function of the curve fl6f). The function F does the job. 

Theorem 2 (Vk \TT\j ) There exists a real function F of the (rescaled) harmonic moments 
t, tk, ik such that 

dF - OF 
CU = vr — , C_ fe = 7T^, k>l. (17) 
dt k dt k 

In particular, this implies the symmetry of derivatives of the inner harmonic moments 
with respect to the outer ones: 

dC_ k = dC^ dCU = dC^ 
dtj dt k ' dtj dt k 



For general Whitham hierarchies, the function F was introduced in ||14|| . In the dis 
persionless Toda case, it is the dispersionless limit of logarithm of the r-function [ITU 



This function obeys the dispersionless limit of the Hirota equation (a leading term of the 
differential Fay identity |T7|, 0): 

{z-0 exp( E — ^ n r m )=^e X pf-E^- fe Vcexp(-E^r fe ] (19) 

This is an infinite set of relations between the second derivatives F nm = d tn d tm F, F 0m = 
d t d tm F. The equations appear while expanding the both sides of (|19D in powers of z and 

C- 

To complete the identification with the objects of the theory of Whitham hierarchies, 
we just mention that the function //(A) (0) is the quasiclasical limit of the Orlov-Shulman 
operator for the 2D Toda lattice. This function obeys the conditions {A,//} = A, 
{A, Jx) = —A. The string equation can be then written as a relation 

A = A"V (20) 

between the Laurent series, together with the reality condition //(A) = /2(A). Written 
in this form, it admits generalizations which are also relevant to the Laplacian growth 
problem, and which we are going to discuss now. 
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Consider the Laplacian growth problem for domains symmetric under the Z n -trans- 
formations z — > e 27Tll / n z, I = 1, ... ,n — 1, where n is a positive integer. In this case 
C m = unless m = (modn), i.e., only the moments C^ n are non-zero. The conformal 
map from the exterior of the unit circle to the exterior of such a domain is given by 
z{w) = (A(w n )) 1//n , where X(w) has the same general structure as in (f|). Equivalently, the 
problem may be posed as the Laplacian growth in the wedge (sector) domain restricted 
by the rays axgz = and argz = 27i/n with the periodic condition on the boundary. In 
the latter case, the conformal map is 

z(w) = (AW) 1 /" . (21) 

The LGE in terms of the conformal map z(w,t) has the same form (|S[) if the time is 
rescaled as t — > t/n. In terms of the A («;,£), one has 

{A, A} = n 2 (AA)^. (22) 
Let us introduce the following notation for the non-zero harmonic moments: 

* = — , «* = % k>l (23) 

rnr irk irk 

(cf. (|XD|) ) . The more general string equation ( p2]) amounts to the relation 

A = n n A-V n (24) 
between the Lax and Orlov-Shulman functions (note that this means \z\ 2 = n/J,(z)). 



This relation is familiar from [10. 12]. It is consistent with the hierarchy and selects 
a solution that, as n > 1, is different from the one selected by (|20"D . The Lax-Sato 
equations ([TT|) - ( |i~4"D for the Lax function A and the Hirota equation (|T9|) hold true as 
they stand provided the times are redefined as in (|23|). The solution to the Laplacian 
growth problem in the channel geometry (in the Hele-Shaw cell) can be obtained from the 
above formulas as a somewhat tricky limit n — ► oo. It is unclear, however, if the existing 
finite-dimensional solutions of the LGE [4-6] will survive after these transformations. We 
should also mention that inclusion of the "no-flux" boundary conditions at the walls of 
the wedge, which state that <91ogz/<91ogu> is real, imposes an extra symmetry in the 
solution and doubles the number of wedges. We will elucidate these questions in the near 
future. 

At last we point out that the above formulas make sense for negative integer n, 
too. In particular, the case n — — 1 describes the internal radial problem when oil is 
inside while water is outside. In this case the internal and external harmonic moments 
are interchanged in their role: the internal moments fll5|) together with the Co become 
independent variables (cf. (|2"3"|)) while the external moments are found as derivatives of 
the prepotential function according to (O). 



These results were reported at the XIII NEEDS Workshop (Crete, June 1999). We 
are very grateful to the organizers for the invitation and for the opportunity to share 
these results in such a fruitful and nice atmosphere. We are indebted to Paul Wiegmann 
for collaboration in [|7|, |Tll and many stimulating discussions. The work of A.Z. was 
partially supported by RFBR grant 98-01-00344. 
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